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I. INTRODUCTION 



Whether we wish to compute the mass spectrum of lattice QCD or the phase structure 
of the extended Hubbard model, we are faced with the same problem — extracting physics 
from a theory to which conventional perturbative methods cannot be applied. To date, the 
most popular approach to these problems has been Monte Carlo evaluation of the Feynman 
path integral. Recently, we introduced an alternative, Hamiltonian-based approach called 
the CORE (Contractor Renormalization Group) approximation OJ and applied it to the case 
of the 1+1- dimensional Ising model. In this paper, we significantly extend the method and 
simplify its implementation. The CORE approach defines a systematic, nonperturbative, 
and computable means of carrying out a Hamiltonian version of the Kadanoff- Wilson |2| 
real-space renormalization group transformation for lattice field theories and lattice spin 
systems. The method relies on contraction and cluster expansion techniques. 

The CORE approximation improves upon other methods of implementing approximate 
real-space renormalization group transformations on Hamiltonian systems in several ways. 
First, our methods make it possible to define a gauge-invariant renormalization group trans- 
formation for any abelian or nonabelian lattice gauge theory, something which was not pos- 
sible in earlier schemes. Second, it is no more difficult to treat fermions than bosons when 
one uses these methods. Third, it is easy to add a chemical potential to the Hamiltonian 
for a system such as the Hubbard model in order to tune the density of the ground state, 
a difficult feat in earlier Hamiltonian real-space renormalization group schemes. Finally, 
CORE allows us to map a theory with one set of degrees of freedom into a theory described 
in terms of a very different set of degrees of freedom but possessing the same low-energy 
physics. Within the context of lattice gauge theories, this means we can start from a theory 
of quarks and gluons and map it into a system in which the effective degrees of freedom have 
the quantum numbers of mesons and baryons. Computing such a transformation within the 
Hamiltonian framework was not possible in earlier methods. 

In addition to the above qualitative improvements, there are also substantial quantitative 
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refinements. For example, earlier attempts to compute the ground-state energy density and 
other properties of the 1 + 1-dimensional Heisenberg antiferromagnet using previous real- 
space renormalization group methods [|J or t-expansion techniques H had difficulty match- 
ing the accuracy of Anderson's || naive spin-wave approximation. We will demonstrate that 
the CORE approximation significantly improves on Anderson's calculation without making 
any large spin approximations. Another example which we discuss is the 1 + 1-dimensional 
Ising model. We will show that an easily implemented CORE computation substantially 
improves upon results from earlier methods. 

We close this section with a brief review of the renormalization group (RG) in order to 
contrast the CORE method to previous RG implementations. Next, in Sec. 0, we state 
without proof the rules for carrying out a CORE calculation. We then illustrate the method 
in Sec. |T I E| by applying the rules in four examples: free scalar field theory with single- 
state truncation, the Heisenberg antiferromagnetic spin chain with two-state truncation, 
the anisotropic Ising model with two-state truncation, and free scalar field theory with an 
infinite-state truncation scheme. The rules are then derived in Sec. [TV} In Sec. [V], two issues 
are discussed: the use of approximate contractors, to establish the connection of the CORE 
approach to earlier methods; and the convergence of the cluster expansion, to demonstrate 
the need for summation via the renormalization group. Finally, future applications to the 
Hubbard and extended Hubbard models and lattice gauge theory with and without fermions 
are discussed in Sec. [V|. We also address the issue of relating the contractor renormalization 
group to the familiar perturbative renormalization group in 4 theory. 

A. Preliminary Remarks 

Physical systems in quantum field theory and statistical mechanics involve a large num- 
ber of degrees of freedom and can usually be described in terms of a local Hamiltonian. 
Conventional wisdom says that when the coherence length of such a system is small, the 
properties of the system depend strongly on the form and strengths of the interactions in the 
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Hamiltonian; whereas, when the coherence length is large, many degrees of freedom behave 
cooperatively and the properties of the system are governed primarily by the nature of this 
cooperation with the detailed form of the Hamiltonian playing only a subsidiary role. 

The renormalization group 0, as formulated by Kadanoff and Wilson, is generally 
thought of as a method for treating systems in which the coherence length encompasses 
many degrees of freedom. This method is based on iteratively thinning the degrees of 
freedom in the problem, an approach which is similar to that followed in hydrodynamics 
wherein the innumerable microscopic degrees of freedom are replaced by a much smaller 
set of spatially-averaged, macroscopic variables, such as the density and pressure. In this 
renormalization group method, the thinning is achieved via a sequence of renormalization 
group transformations. 

While the original formulation of the RG method was done for the partition function or 
its path-integral analogue in field theory, the approach has been extended to Hamiltonian 
systems. The basic idea is to construct a real-space renormalization group transformation, 
r, which maps the Hamiltonian H of a theory defined on some lattice L to a new theory Hi 
defined on a coarser lattice Li in such a way that the new theory has the same low-energy 
physics as the original theory. To extract the low energy physics of the original theory, we 
repeatedly apply the transformation r and generate the sequence of renormalized Hamilto- 
nians: Hi = r(H ), H 2 = t{H\), H 3 = r(H 2 ), . . .. This sequence usually approaches a fixed 
point of r, that is, a Hamiltonian H* satisfying t(H*) = H*. Each renormalized Hamil- 
tonian in this sequence possesses the same low-energy physics, but the degrees of freedom 
have been thinned. Eventually, the number of remaining degrees of freedom lying within 
the coherence length will be small and the resulting Hamiltonian will be more amenable to 
solution. 

Generally, the same transformation r is used for each iteration; however, this is not 
required. The use of different transformations for each iteration is clearly impractical, but 
the use of a different transformation for the first one or few steps could be a powerful 
generalization of the method, facilitating great simplifications. A quantum field theory 
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could be mapped into a generalized spin model; QCD could be mapped into a theory of 
interacting hadrons. 

Defining and carrying out the thinning transformations is the key to the RG approach. 
The RG transformation r is usually defined by requiring invariance of the partition func- 
tion or its path-integral analogue in field theory. The RG method exactly describes the 
low-lying physics as long as r can be exactly implemented, which is rarely the case. In 
practice, approximations must be made, such as those made in the e-expansion f7], the use 
of perturbative matching as in the heavy-quark effective field theory |§ and nonrelativistic 
QCD U, and stochastic estimation as in the Monte Carlo renormalization group [10] ap- 



proach. The CORE approach is a new and powerful method for defining r and computing 
t(H) which relies on contraction and cluster techniques. In contrast to other methods, the 
approximations made in the CORE approach do not limit the usefulness of the method to 
any restricted range of coupling constants or other parameters in the theory. The CORE 
approach works well not only near a critical point when the coherence length is large, but 
also in instances where it is small. It is a general method for solving any lattice Hamiltonian 
problem. 

CORE computations begin by defining the way in which the new lattice is coarser than 
the original lattice. We begin by partitioning the lattice into identical blocks. The Hilbert 
space of states corresponding to each block is then truncated by discarding all but a certain 
number of low-lying states; we generally retain enough states so that the truncated degrees 
of freedom on a block resemble those of a site on the original lattice. The renormalized or 
effective Hamiltonian H ren in this truncated space of states is then defined in terms of the 
original Hamiltonian H by 

H^= lira [T(t) 2 H [T(t)HT(t)j [T{t) 2 ]~i, (1) 

t — >oo 

where T(t) = e~ tH is the contractor and [. . .] refers to truncation to the subspace of retained 
states. There is a one-to-one correspondence between the eigenvalues of the renormalized 
Hamiltonian and the low-lying eigenvalues of the original Hamiltonian. In general, the 



renormalized Hamiltonian cannot be exactly determined; CORE approximates H mn using 
a finite cluster expansion, an approximation which can be systematically improved. Matrix 
elements of various operators can also be evaluated in CORE by defining a sequence of 
renormalized operators. 

We use the phrase "CORE technology" to refer to the set of tools which allow us to 
systematically and nonperturbatively compute an arbitrarily-accurate approximation to the 
exact renormalization group transformation for a lattice field theory or spin system without 
having to diagonalize the original infinite-volume theory. The power of these methods is 
that usually only a few terms in the cluster expansion of the renormalized Hamiltonian yield 
remarkably good results. 



In this section, we state, without proof, the rules for carrying out a CORE computation. 
We assume that we are studying a theory defined by a local Hamiltonian H on a regular 
lattice of infinite extent in some number of dimensions. 

A CORE computation proceeds as follows: 

1. First, divide the lattice into identical, disjoint blocks Bj. Denote the space of states 
associated with block Bj by TLj and denote the common dimension of each of these 
spaces by N. 

2. Define a truncation scheme by selecting a low-lying subspace Vj C TLj of dimension 
M < N on every block; the same subspace should be chosen on each block. In what 
follows, we will denote the M retained states by \(/> s (j)) and use them to construct the 
projection operators 



II. THE RULES 
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Let [. . .] denote truncation to the subspace spanned by the taking tensor products 
of the states \<f> s (J))- Thus, for any operator O, the truncated operator is defined 
as [O] = P^OP. Note, choosing to retain states such that the truncated degrees of 
freedom on a single block resemble those associated with a single site on the original 
lattice ensures that the renormalized Hamiltonian will take a form similar to that of 
the original Hamiltonian, facilitating the iteration process; however, sometimes it is 
useful to make a different choice and map the original theory into one formulated in 
terms of new degrees of freedom. 

3. Compute (see below) the renormalized Hamiltonian defined in Eq. [I], H[ en = t(H), 
and the renormalized operators corresponding to any matrix elements of interest. 

4. Repeat the above steps using H™ n to obtain = r(H™ n ). Iterate this process 
until the renormalized Hamiltonian is simple enough that its low-lying eigenvalues can 
be determined. 



Because the Hamiltonian is extensive (a concept we will define later) and the block-by- 
block truncation preserves this property, the renormalized Hamiltonian can be approximated 
using the finite cluster method (FCM). This method was first used by Domb fll| in the 
application of the Mayer cluster integral theory to the Ising model. A formal proof of 
the method in the Ising and Heisenberg models was then presented by Rushbrooke |TJ| . 



The method was later generalized by Sykes et al. [13]. The finite cluster method expresses 
any extensive quantity in an infinite volume as a sum of finite-volume contributions. The 
procedure is simple to implement and provides numerous means of detecting computational 



errors. A general statement of the method can be found in Ref. ||T4| 



Evaluation of H ren by the finite cluster method is accomplished in the following sequence 
of steps: 



1. Compute the renormalized Hamiltonian for a theory defined on a sublattice which 
contains only a single block Bj (how this is done will be described below). Denote this 
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Hamiltonian by H r (Bj) = hi(Bj). This yields all of the so-called range-1 terms in the 
cluster expansion of the renormalized Hamiltonian. 

2. Calculate the renormalized Hamiltonian H r (Bj, Bj + i) for a theory defined on a sub- 
lattice made up of two adjacent {connected) blocks Bj and Bj + \. The range-2 con- 
tributions to the cluster expansion of the renormalized Hamiltonian on the infinite 
lattice are obtained by removing from H T (Bj, Bj +1 ) those contributions which arise 
from terms already included in the single block calculation: 

h 2 (B j: B j+1 ) = H r (B j: B j+1 ) - h^Bj) - (4) 

3. Repeat this procedure for sublattices containing successively more connected blocks. 
For example, for a sublattice consisting of three adjacent blocks Bj, Bj + i, and Bj +2 , 
use 

h 3 (Bj, B j+1 , B j+2 ) = H r (Bj, B j+1 , B j+2 ) - h 2 (Bj, B j+1 ) - h 2 (B j+1 , B j+2 ) 

-h^-hiBj+J-hiBw). (5) 

Since the renormalized Hamiltonian is extensive, only connected sublattices need to 
be considered. Recall that a quantity is extensive if, when evaluated on a discon- 
nected sublattice, it is the sum of that quantity evaluated separately on the connected 
components of the sublattice. A truncated cluster expansion can then be defined by 
neglecting clusters larger than some specified range. 

4. To complete the determination of the renormalized Hamiltonian H ren on the infinite 
lattice, sum the connected contributions h r (Bj, . . .) from the finite sublattices accord- 
ing to their embeddings in the full lattice. For example, on an infinite one-dimensional 
lattice: 

oo oo 

# rcn = E EM/;, n,-r ,!. (6) 

j=— oo r=l 
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Express H ren in terms of block- variables such that the form of H ren resembles that of 
the previous Hamiltonian in the RG sequence. 

The key ingredient of CORE is the method used to explicitly construct the renormalized 
Hamiltonian H r (G) and other renormalized operators on a given cluster or sublattice G. 
While, in principle, the appropriate generalization of Eq. [1] completely specifies what has 
to be done, in practice, an attempt to compute this quantity by brute force will run into 
problems since the operator [[T(i)]] becomes singular as t — > oo. To see how this problem 
arises, consider a sublattice comprised of R connected blocks B\. . . Br. Let H(G) denote the 
Hamiltonian obtained by restricting the infinite lattice H to the sublattice G and suppose 
that we truncate to the subspace V — V\ ® ■ ■ • ® Vr spanned by the M R states \<f> a (G)). 
Remember that the states \<f> a (G)) are tensor products of the retained states on each of the 
R blocks in the cluster G. Let us denote by \ep(G)) the eigenstates of H{G) with eigenvalues 
€p(G) and expand the states \<p a (G)) in terms of these eigenstates: i.e., 

\UG)) = H^MG)). (7) 

It then follows that 

T(t)\MG)) =Y,^e-^\ep(G)), (8) 

P 

from which we see that all states \<p a (G)) which have a nonvanishing overlap with the ground 
state of H(G) contract onto the same state as t — > oo. This causes great difficulties if we 
attempt to numerically compute [[T(t) 2 ]] -1 / 2 . Fortunately, there is an elegant and simple 
solution to this problem which avoids explicit computation of [[T(t) 2 ]] -1 / 2 : make a unitary 
(or orthogonal) change of basis, S(G), on the states \4> a (G)) such that each state in the 
new basis contracts onto a unique eigenstate of H(G). In this new basis, the computation of 
H r (G) is then straightforward. The discussion which follows specifies the rules for computing 
the necessary change of basis and H r (G) for a general H(G). We state these rules in full 
generality so as to allow for the special situation in which H (G) has degenerate eigenvalues, 
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and then apply them to successively more complicated examples in order to show how they 
work in practice. 

H I (G) and the change of basis S(G) may be determined as follows: 

1. Find the eigenstates \ep(G)) and corresponding eigenvalues ep(G) of H(G), where 
(3 — . . . N R — 1. Order these states so that ep < ep + \. 

2. Construct the M R x N R matrix 

Q{G) aP ={<f> a (G)MG)). (9) 

Each row of Q(G) gives the expansion of one of the retained states in terms of the 
eigenstates of H{G). Each column of Q(G) gives the projection of some eigenstate 
into the truncated subspace. Also, let S(G) be the M R x M R identity matrix. Set 
m = M R , p = 0, and q = 0. 

3. Copy the first g columns of Q(G) into an m x g matrix C, where g is the degen- 
eracy of the lowest-lying eigenvalue. If the ground state of the cluster is nonde- 
generate, then g = 1. The columns of C correspond to the degenerate eigenstates 
\e q (G)), ■ ■ ■ , |e g+g _i(G)). Having formed C, perform a singular value decomposition 
(SVD), writing 

C = UZV\ (10) 

where U is an m x m unitary matrix, V is a g x g unitary matrix, and S is an m x g 
matrix of the form 



(11) 

A = diag(cri,---,cr r ), (12) 



i 0( m _ r )xr 0(m—r)x(g— r) j 



where the elements cij are real and satisfy <J\ > <J2 > • • • > cr r > and r < min(m, g) is 
the rank of the matrix C. In other words, use the SVD to construct orthonormal 



10 



bases for the nullspace and range of the matrix C. Note that the SVD theorem 
guarantees that such a decomposition exists and that S is unique. 

4. Multiply U* Q{G) and, by abuse of notation, once again call the result Q(G). Then 
discard the first g columns and the first r rows of the new Q(G). The resulting matrix, 
which we again call Q(G), is now an (m — r) x (N R — q — g) matrix. Note, r may be 
zero. 

5. Form the matrix 

d o \ 



R 



L pxp "pxm 
y u mxp u mxm J 



(13) 



and multiply i? S(G); call the result S(G). Define the states 

IWi(G)> = E ^ ( 14 ) 



s'=l 



with corresponding degenerate eigenvalues 7^ +s _i(G) = e q (G), for s = l...r. Set 
p — > p + r, q ^ q + g, and m — > m — r. 

6. Repeat steps ^| to |5| with higher and higher energy eigenvalues until p = M R . In step 
^, g is now the degeneracy of the lowest-lying remaining eigenvalue. At the end of this 
procedure, we will have constructed a unitary M R x M R matrix S(G) and a set of 
eigenstates {\1p(G))} with energy eigenvalues Tp{G) for (3 = 1 . . . M fi . 

In the discussion which follows, it will be convenient to make the following definitions: 

Definition: The eigenstates \Tp(G)) are referred to as the remnant eigen- 
states of H(G) in V. The set of these M R remnant eigenstates is called the 
contraction remnant. The matrix S(G) is referred to as the triangulation 
matrix. 



As we already noted, the triangulation matrix S(G) is simply a change of basis, taking 
us from the original basis {\4> a (G))} of retained tensor-product states in the truncated 

11 



subspace to another basis {|£a(Cr))} in which only the first state has a nonvanishing 
overlap with the ground eigenstate, only the first and second states have nonzero 
overlaps with the first excited eigenstate, and so on; hence, S(G) a p = (£ a (G)\(j)p(G)). 
The remnant eigenstates are essentially the M R lowest-lying eigenstates of H{G) whose 
projections into V are nonvanishing and cannot be written as linear combinations 
of lower-energy eigenstates projected into V. In other words, the projections of the 
remnant eigenstates into V are all linearly independent. Within degeneracy subspaces, 
the eigenstates must be rotated in order to eliminate all linear combinations whose 
projections in V are zero or completely expressible in terms of the projections of lower- 



lying eigenstates. Note that the singular value decomposition theorem |15[ guarantees 
the existence of the triangulation matrix and the contraction remnant. 

7. In the basis of the remnant eigenstates, construct the matrices 

Ht{G)c0 = (T a (G)\H(G)\Tp(G)) = 6 aP T a (G), (15) 
T (GU = (T a {G)\0{G)\T p {G)), (16) 

where 0(G) is some operator of interest defined on the sublattice G. 

8. The renormalized operators given in terms of the triangulation matrix and 
the operators evaluated in the contraction remnant: 

H V {G) = S^{G) H r (G) S{G), (17) 
T {G) = S^{G) O r {G) S{G). (18) 



Note that the CORE approach described here differs from that described previously [|TJ . 
In our earlier formulation of the method, the contractor T(t) in Eq. [I] was approximated by a 
product of exactly computable exponentials. The variable t was then treated as a variational 
parameter, adjusted so as to minimize the mean-field energy in each RG iteration. 
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III. FOUR EXAMPLES 



To better illustrate the method and demonstrate its effectiveness, we now apply these 
rules in four examples. Each of these examples, free scalar field theory with single-state 
truncation, the 1 + 1-dimensional Heisenberg antiferromagnet with two-state truncation, 
the 1 + 1-dimensional Ising model with two-state truncation, and free scalar field theory 
with infinite-state truncation, has been chosen to clarify a particular aspect of the rules. 



A. Single-State Truncation: Free Scalar Field Theory 

First let us discuss a massless (// = 0) free-field theory. Free scalar field theory on a 
lattice is just a set of coupled harmonic oscillators, 



VjT + ^0(jT + + 1) - <Kj)) 2 



(19) 



2 w 2 r w, 2 

where [<f>(j),H(k)] = i5jk- The simplest possible truncation procedure we can adopt is to 
keep the number of sites fixed and truncate to a single state per site. Begin by dividing H 
as follows: 

H = J2(H(J) + V(j)), (20) 

3 

H(j) = \ (n(j) 2 + 20(j) 2 ) , (2i) 

nj) = -0(j)0(j + l)- (22) 

Truncate by keeping only the ground-state of H(j) for each site j; i.e., keeping the oscillator 
state \uj(j)) of frequency uo = \/2. Note, this procedure truncates the entire Hilbert space 
to a single product-state and therefore the renormalized Hamiltonian will be a 1 x 1-matrix, 
as will each term in the expansion 

H Icn = 5> r (j) conn . (23) 

3,r 

Since the CORE procedure guarantees that H ren has the same low energy structure as the 
original theory, keeping only one state means that we will only be able to compute the 
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ground-state energy of the free scalar-field theory. We will see that all of the terms h r (j) conn 
are independent of j and so it follows from Eq. |23| that the ground-state energy density will 
be given by 

oo 

^ee-ficld = E^O') COnn (24) 
r=l 

for any fixed j. 

Following the basic rules, truncate H(j) to obtain 

MjT onn = \V2, (25) 

where hi(j) conn can be thought of as either a 1 x 1-matrix or as a c-number. 

To compute the range-2 contribution to the energy density, we must diagonalize the 
two-site Hamiltonian 

# (jVsite = \ (n(j) 2 + 20(j) 2 + n(j + 1) 2 + 200' + 1) 2 ) - <f>U)<f>U + 1), (26) 

and expand the tensor product state \oj{j)) (g> + 1)) in terms of the eigenstates of 
if (j) 2 -sites- Since this tensor-product state has the exact two-site ground state appearing in 

(2) 

its expansion in terms of the two-site eigenstates, H T (j) is a 1 x 1 matrix whose single entry 
is the exact ground-state energy of i?2-sit e ; E 2 = |(a/3 + 1). Furthermore, since S has 
to be a 1 x 1 orthogonal matrix it is trivial. It follows from these facts that the connected 
range-2 contribution to the ground-state energy density is given by 

h 2 ( 3 y o ™ = E 2 (j) - 2Mj) c ° nn = ^(V3 + 1 - 2V2). (27) 

To construct the range-3 term, find the ground-state energy of the three-site problem, E%, 
and then subtract twice the range-2 contribution, because we can embed a connected two- 
site sublattice in the three-site lattice in two ways, and three times the range-1 contribution, 
because the single-site can be embedded in the three-site sublattice in three ways; i.e., 

h c r n U) = EsU) - 2/i 2 (j) c ° nn - 3/iiOT olin (28) 
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To compute the range-r contributions, find the exact ground-state energy of the r-site Hamil- 
tonian, E r , and then subtract the lower order s-range connected contributions as many times 
as the corresponding connected s-site sublattice can be embedded in the r-site problem: 

h r (jr™ = E r (j) - 2/ V _ 1 (j) conn - 3/v_ 2 (j) c ° nn - ... - rh^Y^. (29) 

We wish to emphasize the unusual nature of this formula in that we calculate the energy 
density of the infinite- volume Hamiltonian system by exactly solving a series of finite-lattice 
problems, each defined with open boundary conditions, and recombine these results to cancel 
out finite-volume effects. The results shown in Table | show the way in which the partial 
sums 

n 

e n = Y,hrti) comi (30) 

r=l 

converge to the true ground-state energy density. The surprising result, given that the 
energies E r are computed for problems with open boundary conditions, is that the finite- 
volume effects appear to cancel to order 0(l/r 3 ), rather than 0(l/r) as one would expect 
for a theory defined on a finite lattice with open boundary conditions, or like 0(l/r 2 ) which 
one would expect for a theory defined with periodic boundary conditions. At this time we 
do not completely understand why the convergence is this rapid, but this behavior is seen 
in all of the examples we have studied. 



B. Two-State Truncation: Heisenberg Antiferromagnet 

There are several reasons for studying the Heisenberg antiferromagnet. First, the model 
exhibits spontaneous breaking of a continuous symmetry in two and three spatial dimensions 
and, although in one spatial dimension the Mermin- Wagner [|16|] theorem forbids a nonva- 
nishing order parameter, the theory still has a massless particle; it is interesting to see if 
we can obtain the ground-state energy density, the massless spectrum, and the vanishing of 
the staggered magnetization by means of a simple CORE computation. Second, this theory 
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is exactly solvable by means of the Bethe ansatz [|T7| and so we can compare our results 
to the exact ground-state energy density e exa ct = — l n (2) + 1/4 = —0.443147. Third, there 
is a computation by Anderson, based on an approximate spin-wave computation, which 
reproduces the spin-1/2 antiferromagnet energy-density to within 2.5%. Although this ap- 
proximate result is based on treating the spin-1/2 system as if it had spin-iV, for N ^> 1, and 
then evaluating the result for N = 1/2, it has been difficult to do as well by earlier Hamil- 
tonian real-space renormalization group methods; we are finally able to exhibit a simple 
approximate CORE computation which does significantly better than Anderson's spin-wave 
computation working with spin-1/2 from the outset. The final reason for studying this case 
is that the symmetry of the model makes it possible to describe the details of the compu- 
tation in a straightforward manner. In particular, it is simple to explain the need for, and 
construction of, the triangulation transformation S which we referred to when we stated the 
basic rules for doing a CORE computation. 

The Heisenberg antiferromagnet is a theory with a spin-1/2 degree of freedom s(j) at- 
tached to each site j of a one- dimensional spatial lattice and a nearest-neighbor Hamiltonian 
of the form 

H = J2s(j)-s(j + l). (31) 

3 

The s(j)'s are operators which act in the single-site Hilbert spaces Tij and satisfy the familiar 
angular momentum commutation relations 

Mi), sp(l)] = iSjiea^s^j). (32) 

To analyze this problem, divide the lattice into three-site blocks and label each block by 
an integer j. The sites within each block are labelled by the integers {3j, 3j + 1, 3j + 2}. 
Corresponding to this decomposition of the lattice into blocks, divide the Hamiltonian into 
two parts, H B and Vbb, 

#b = E MJ) = E [ s ~T 3 j) ■ S ~T 3 J + !) + S ~*( 3 J + !) • S ~*( 3 J + 2 ) ] : ( 33 ) 

3 3 

Vbb = E ^bb(j) = E s "(3j + 2) • s(3j + 3). (34) 

3 3 
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Truncate by keeping the two lowest-lying eigenstates of H&(j) for each block B(j) so as to 
produce a new coarser lattice which again has a spin-1/2 degree of freedom associated with 
each of its sites. Diagonalizing H#(j) is a simple exercise in coupling three spins; i.e., 

H B (j) = s(3j) ■ *(3j + 1) + s(3j + 1) • s(3j + 2), (35) 
= s(3j + 1) ■ (s(3j) + s(3j + 2)) , (36) 
= l{sL(j)-Sf 0+2) (j)-3/A) J (37) 

where S tot (j) = s{3j) + s{3j + 1) + s(3j + 2) and S 0+2 (j) = s{3j) + s(3j + 2) . From Eq. |7| we 
see that the eigenstates of -^b(j) can be labelled by the eigenvalues of S? ot (j) and S^ 0+2 ^(j), 
and the two lowest-lying eigenstates belong to the spin-1/2 multiplet for which the spins on 
sites 3j and 3j + 2 couple to spin-1. We denote these two degenerate states by | fj) and 
| lj) and use them to construct the projection operator 

p = U p U); PU) = \WC\i\ + \ii)(iil ( 38 ) 

Using the P(j)'s we construct the connected range-1 operators 

Mj) con„ = l> U) lf ]>>U) l> U) = (39) 

where l(j) stands for the 2x2 identity matrix. 

To obtain the connected range-2 term ^2(j) com \ construct the Hamiltonian for the two- 
block or six-site problem. Since this Hamiltonian commutes with the total-spin operators 
for the six-site sublattice, the eigenstates of H six _ sites will fall into spin-3, spin-2, spin-1 or 
spin-0 multiplets. The following state, 

^ (I tili+i) - I liti+i)) , (40) 

is the unique linear combination of the original tensor-product states which has total spin 
zero; hence, only spin-0 states appear in the expansion of this state in terms of eigenstates 
of if S ix-sites- The lowest-lying eigenstate of if S i x _ S it es appearing in the expansion of this spin- 
state is the ground-state of H six _ sitcs whose eigenvalue we denote by e . Similarly, the 
following states 

17 



tjtj+i) > 7/| (I Tiij+i) + Ujti+i)) > I 



(41) 



are linear combinations of the original tensor-product states which have total spin 1 and total 
z-component of spin M z = +1,0,-1, respectively. The lowest-lying eigenstate of if S i x _ S ite S 
appearing in each of these spin-1 combinations is that member of the lowest-lying spin-1 
multiplet having the appropriate value of M z ; hence, each of these states contracts onto a 
unique eigenstate of -ff S i x _ S ites- If we denote the degenerate eigenvalue of these eigenstates 
by ei, then the operator H^\j) has the form 



(42) 



( e ^ 
ei 
e 1 
ei 

using these remnant eigenstates as our new basis states. We could use the explicit form 
of the triangulation matrix S, which rewrites the original tensor product states in terms of 
these spin eigenstates, to transform this back into the original tensor-product basis 
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tjtj+i); I tj!j+i)> I Ijtj+i)) I ijlj+i); 



(43) 



but this is unnecessary since symmetry considerations require S^H^\j)S to have the form 



(44) 



Eq. ^ can be rewritten in terms of the total spin operator for sites j and j + 1 to obtain 

(45) 



3 1 

6 = /3o-tA ; ei = A) + -/%.. 



While the symmetry of this system makes it possible to determine e and e x analytically, 
it is more convenient to compute it numerically. To six significant figures, this calculation 
yields 



e = -2.493577 ; e x = -2.001995. 



(46) 



To construct the connected range-2 term, we subtract the two ways of embedding the 
one-block sublattices into the connected two-block sublattice 

h 2 (jy oim = (A) + 2)1, ® + Mi) ■ 8(3 + 1). (47) 

We could go on to compute range-r connected terms for r > 2, but we will stop at range-2 
and define the approximate renormalized Hamiltonian by 

H rcn = ]T Clj + p 1 s(j) ■ s(j + 1), (48) 

j 

where C = /3 + 1. Clearly this approximate Hamiltonian, except for the trivial addition of 
a multiple of the unit matrix, has the same form as the original Hamiltonian. When this 
happens, we say that the theory is at a critical point, and \j3i\ < 1 implies that it has no 
mass gap. ( The logic which says that \j3i\ < 1 implies no mass gap is that if we iterate 
the renormalization group transformation, then eventually only the c-number part of the 
Hamiltonian will remain. Since the interaction part becomes vanishingly small, eventually 
all of the low-energy states of the theory must have a vanishingly small energy splitting. 
Hence, the theory must have a vanishingly small mass gap.) 

To extract the ground-state energy density, we have to pay attention to the constant 
term. After the first transformation, we see that this term will make a contribution to the 
ground-state energy density equal to C/3, where the factor of 1/3 appears because each 
site on the new lattice corresponds to three sites of the old lattice. Remembering this and 
performing the renormalization group transformation on the interaction term (3is(j) -s(j + l), 
we generate a new renormalized Hamiltonian of the form 

HL = E PiClU) + (%3U) ■ s(j + 1). (49) 

j 

Accumulating the new constant f3\ C/9 into the previous computation of the energy-density 
(where the 1/9 comes from the fact that one point on the new lattice corresponds to nine 
points on the original lattice), we again have a new Hamiltonian which has the same form 
as the original Hamiltonian, except that it is multiplied by the factor f3\. Repeating this 
process an infinite number of times yields a series for the ground-state energy density 
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n 



c 



£ 



Ten— grp 



3(1 -A/3) 



-0.4484462, 



(50) 



which agrees well with the exact result £ e xact = —0.443147. Thus, this simple range-2 
calculation gives a result which is good to about one-percent; this is more than a factor of 
two better than that obtained from Anderson's spin-wave computation. Note that this very 
simple calculation yields the exact mass gap. One also finds that the staggered magnetization 
M. = Z)j(— l) ,J 's«(j) vanishes (note that one obtains a nonvanishing staggered magnetization 
on two- and three-dimensional spatial lattices). 

This completes our present discussion of the antiferromagnet. We will return to it again 
in the section on questions of convergence since it has something to teach us about the 
reliability of single-state truncation calculations. 



We now revisit the 1 + 1-dimensional Ising model which we discussed in Ref. JIJ using 
an earlier formulation of the CORE approximation. While our earlier treatment was quite 
successful in extracting the physics of the model, our new approach produces better results, 
is less computationally intensive, and is much easier to implement and explain. There are 
two main reasons for treating this example in some detail. First, remarkably accurate results 
can be obtained even when considering only terms up to range-3 in the renormalized Hamil- 
tonian. Secondly, this problem does not have the high degree of symmetry of the Heisenberg 
antiferromagnet and so the construction of the operator S must be done explicitly. 

The Hamiltonian of the 1+ 1-dimensional Ising model is 



where j labels the sites on the infinite one-dimensional spatial lattice and < A < 1. This 
model is interesting for several reasons. First, it exhibits a second-order phase transition 



C. Two-State Truncation: The 1 + 1 dimensional Ising Model 



= - [CA^O') + s A^(j>x(j + 1)] , 



(51) 



j 

c\ = cos(A7r/2), s\ = sin(A7r/2), 
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at A = 1/2; for A < 1/2, the ground state of the system is unique, the order parameter 
(<r x (j)) vanishes and the excited states are localized spin excitations; when A > 1/2, the 
ground state is twofold-degenerate corresponding to values of the order parameter (a x (j)) = 
±(1 — cot 2 (A7r/2)) 1//8 and the excitations are solitons (kinks and antikinks). Secondly, the 
model is exactly solvable and so we have exact results with which to compare. Thirdly, 
the model has much less symmetry than the Heisenberg model and so the structure of the 
renormalization group transformation is richer. 

In order to show how a more complicated approximate renormalization group transfor- 
mation works, we once again adopt a two-state, three-site block truncation algorithm, but 
now we compute h r (j) conn for r = 1, 2, 3. Because the Hamiltonian is more complicated than 
that of the antiferromagnet and computing the connected range-3 terms involves solving for 
the eigenstates of a nine-site problem, we must resort to numerical methods to carry out 
the computation. We numerically diagonalize the 512 x 512 nine-site Hamiltonian matrix. 
However, since we only need a few low-lying states to compute S, we could significantly 
reduce the computational cost by using the Lanczos method. While unnecessary for this 
simple problem, the application of the Lanczos method to the construction of S will be very 
useful when studying more complicated theories. 

The starting Hamiltonian is invariant under parity and the simultaneous transformation 
s x(j) —> — s x(j)- Our thinning algorithm preserves this symmetry so that the most general 
form the renormalized Hamiltonian can take is 

H™ = -^c a O a (j), (52) 
O a {j) = cr ao (j) a ai (j + 1) • • • a ar (j + r), (53) 

where the C a S £1X6 the couplings, a labels the different types of operators which can appear, 
and j is a site label. Given the symmetries of the original Hamiltonian which will be pre- 
served in the renormalized Hamiltonian, we see that the only two possible one-site operators 
are = {u,z}, where u denotes the identity operator; in other words, the only one-site 
operators are O u (i) = a u {i) = 1 and O z (i) = a z (i). Similarly, the only two-site operators 
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which are consistent with the symmetries of the problem are = {xx,yy, zz}, and the 
only three-site operators which can appear are a^ 3 ^ = {xzx,xux,xxz, zxx,yzy,yuy,yyz, 
zyy, zuz, zzz}. 

Since the original form of the Hamiltonian given in Eq. [51] is just a special form of 
Eq. we will discuss the truncation procedure for the general case. Once again we work 
with blocks B,, containing the points {3j, 3j + 1, 3j + 2} and keep the lowest two eigenstates 
of the generic block Hamiltonian 

H B (j) = -c z [O z (3j) + O z (3j + 1) + O z (3j + 2)] - c xx [O xx (3j) + O xx (3j + 1)] 

- c yy [O yy {3]) + O yy {3j + 1)] - c zz [O zz {3]) + O zz (3j + 1)] - c xux O xux {3j) 

CyuyOyuy{3j^) C zuz Q zuz ( 3 J ') C xzx (D xzx (3j^ C xxz (D xxz (3j^) C zxx O zxx ( 3 J ') 
CyzyOyzy(3j^) CyyzOyyz(3j^) CzyyO zyy (3j ') C zzz O zzz (3j^ . (^4) 

When we truncate H-q{j) to these two states, we obtain the new range-1 terms. Since 
hi(j) conn is a diagonal matrix in this basis, it can be written as a sum of a multiple of the 
unit matrix and a multiple of s z (j). 

If we denote the eigenstates of the three-site block Bj by | 1j) and | jj), then the con- 
nected range-2 contributions to the renormalized Hamiltonian are obtained by first solving 
the two-block or connected six-site problem exactly and expanding the four tensor-product 
states ^(B^Bj'+i)) = { | fjfj+i)> | fj!j+i), I ljfj+i) and | in terms of the 

eigenstates \e^(Bj, B J+ i)) of this problem. To compute S for the range-2 problem, begin 
by constructing the 4 x 64 overlap matrix Q a p = (<f) a (Bj, Bj + i)|e 1 a(Bj, B J+ i)). Note that 
each row of Q gives the expansion of each of the four tensor-product states in terms of 
the eigenstates of -fZ" six _ si t e . Ensure that the eigenstates are arranged in order of increasing 
eigenenergy. 

The construction of S now proceeds iteratively. Begin by focusing attention on the first 
column of Q; this is a 4 x 1 matrix C\ whose entries contain the overlaps of the four tensor- 
product states with the eigenstate of lowest energy. If C\ has any nonzero entries, then 
we can find a rotation matrix R\ such that C\ can be brought into a form where only its 
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upper entry is nonzero. Finding such an R\ is equivalent to constructing the singular value 
decomposition of C\. Using R ly transform Q to Qi = R 1 Q and then focus attention on 
the 3 x 63 submatrix obtained by eliminating the first row and first column of Q\\ call the 
resulting matrix Q\. 

Apply the same reasoning to Q\. Focus on the first column of Qi, denoted by C 2 . If C 2 
contains some nonvanishing entries, construct the orthogonal 3x3 transformation U 2 which 
brings C 2 into the standard form where only the upper element is nonvanishing. Again, this 
is equivalent to performing the singular value decomposition of C 2 . Next, define a 4 x 4 
matrix R 2 as 



R? 



^io o o x ' 



ul 



(55) 



/ 



Transform Q to Q 2 = R2R1Q, then construct the 2 x 62 submatrix Q 2 formed by eliminating 
the first two rows and columns from Q 2 . 

Next, construct the 2x2 matrix U\ which brings the first column of Q 2 into standard 
form, extend U\ to a 4 x 4 matrix 



i* 3 



V 



(56) 



^100 o s 

10 

00 

00 ; ul 

and then define the triangulation matrix = R^R 2 Ri- Also, define the diagonal 4x4 

(2) 

Hamiltonian H)-' from the four lowest eigenvalues of H six _ sitcs . Then the 4x4 connected 
range-2 operator in the renormalized Hamiltonian is given by 

h 2 (j) conn = S m H ( ^S (2) - /ii(j) c ° nn - ^{j + l) conn . (57) 

Note that we have simplified this discussion by assuming, as is usually the case, that 
the eigenvalues of if S i x - S ites are nondegenerate and that the four tensor-product states have 
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nonvanishing overlaps with the four lowest-lying eigenstates. If this is not the case, then we 
have to generalize this discussion slightly. In the event that some eigenstates do not occur 
in the expansion of the tensor product states, the corresponding matrix Q\ or Q 2 will have a 
first column in which all entries are zero. When this happens, simply eliminate this column 

and use the first nonvanishing column to define the rotation matrix; the corresponding 

(2) 

eigenvalue is then used in Hj- . When an eigenvalue is (7-fold degenerate, include in C\ or 
C*2 all g columns of Qi or Q 2 corresponding to the eigenvectors in the degeneracy subspace 
and then carry out the singular value decomposition Cj = UjT.jVj . The required rotation 
matrix is again obtained from Uj, but now Vj is needed to construct the remnant eigenstates 
from the degeneracy subspace. Taking degeneracies into account can become important after 
a large number of iterations as the renormalized Hamiltonian flows closer and closer to one 
of its fixed points. 

Constructing the range-3 connected terms proceeds the same way, except now we have 
to work with three adjacent blocks Bj, B J+1 and Bj +2 . Now the matrix Q has eight 
columns corresponding to the tensor-product states | Tjtj+itj+2)) I ljtj+itj+2), I Tj-lj+iTj+2) ? 
I hh+dj+2), I lilj+iTj+2), I liTi+ilj+2), I hij+iij+2), and I ijij+iij+2), and we must 
compute (assuming nondegeneracy of the spectrum and no missed states) the eight lowest 
eigenstates of -ffninc-sitcs(j) in order to construct and Except that we are dealing 
with slightly larger matrices, we go through the same steps described above. The range-3 
contribution to the renormalized Hamiltonian is then given by 

h 3 {j) conn = S^ 3) H^\j)S {3) - /ii(j) coim - h t {j + l) conn - h x {j + 2) conn 

- /i 2 (j) c ° nn - h(J + l) conn - (58) 
Given h^(j), h^ 2 \j), and h^(j), the approximate renormalized Hamiltonian on the lattice 
with one-third as many sites as the original is 

H^ = J2[hi(j) + h 2 (j) + h 3 (j)]. (59) 
j 

Our results were obtained by choosing a specific value of A in the special form of the 
Hamiltonian given by Eq. |5l| in which only and c^j. differ from zero. We then apply 
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the above range-3 CORE procedure. After the first RG transformation, we obtain a new 
Hamiltonian comprised of all allowed operators with nonvanishing couplings; however, most 
of the couplings are small. We iterate the process, obtaining a sequence of renormalized 
Hamiltonians in which the couplings flow until finally, all but one of the coefficients vanish; 
i.e., until one reaches a solvable fixed-point Hamiltonian. Our numerical computations 
show that there are only two possible fixed-point Hamiltonians: one in which only 4°°) is 
nonvanishing, and one in which only is different from zero. 

In Fig. H|, we plot the fractional error in the CORE estimates of the ground-state energy 
density. The dotted curve shows the results obtained in Ref. |3J using the earlier version of 
the CORE approximation. The critical value A c separating the spontaneously-broken phase 
from the unbroken phase is found to be A c ~ 0.50365 which agrees well with the exact value 
of 1/2. 

Extracting the mass gap as a function of A is easily done since both fixed-point Hamilto- 
nians are exactly solvable. Below the phase transition where 4°°) is the only nonvanishing 
coefficient, eigenstates of the Hamiltonian are tensor products of eigenstates of cr z {j), and 
so the mass gap is equal to 24°°' ) - Above the phase transition, the only nonvanishing coeffi- 
cient is 4z and so eigenstates of the Hamiltonian are products of eigenstates of cr x (j). In 
this case, there are two degenerate ground states; the discrete symmetry a x (j) — > —cr x {j) 
is spontaneously broken. In this phase, the low-lying eigenstates are kinks which have mass 
2c < ££'. The results of the CORE computations for the mass gap are shown in Fig. |2]. 

Finally, the magnetization was also studied. A sequence of renormalized magnetization 
operators was computed along with the renormalized Hamiltonian; the starting operator in 
this CORE sequence was M = J2j a x{j)- The renormalized magnetization M ren has a cluster 
expansion given by 

00 

M reil = EE^0T onn , (so) 

j r =l 

where the connected range-r operators m r (j) conn are computed from the truncated one, two, 
and three-block operators, 
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mi(j) co„n = M (l) (j); (61) 

m 2 (j) conn = S m M^\j)S {2) - m^j) 001111 - m x {j + l) conn , (62) 
m 3 (j) conn = S m M^\j)S {3) - m 2 (j) conn - m 2 (j + l) conn 

- rmUr^ - mi (j + l) conn - mi (j + 2) conn , (63) 

where Af!f\j) consists of the matrix elements of in the basis of remnant eigenstates, 
and is the restriction of the full magnetization operator to the r-block sublattice. A 
comparison of the CORE estimates of the magnetization with the exactly known results is 
shown in Fig. |3|. 

Two procedures were used to extract the critical exponent from these calculations. Both 
procedures attempt to fit the logarithm of the magnetization to the form of the exact answer, 
namely: 

ln(M(A)) 

^ ' ~ ln(l - tan(A c 7r/2) V tan(A7r/2) 2 ) ' ( ' 

where ln(M(A)) stands for the logarithm of the computed values of the magnetization and 
A c stands for that value of A at which the theory changes phase. If we attempt to extract 
p(X) by fixing A c = 0.50365, the value above which the CORE computation changes from 
having ^ to c^' ^ 0, then the values of p(X) obtained from this procedure do not lie 
on straight line. Moreover, the average value of p(X) lies between 0.10 — 0.11, which is not a 
very good fit to the exact value 1/8. If, on the other hand, we vary A c and determine its best 
value by fitting the resulting values for p(X) to a straight line, then we obtain a very good fit 
for A c ~ 0.498 and find that p(X) lies in the range 0.1236 < p(X) < 0.126. The discrepancy 
between A c = 0.498 and A c = 0.50365 gives a priori evidence, without knowledge of the 
exact solution, that the determination of the critical point must have an error of about one 
percent due to an accumulation of numerical errors and limiting the computation to range-3 
terms. Fig. [| displays three plots of p(X) for A c = 0.496, 0.498, and 0.500; the best fit to a 
straight line is given by the middle curve which corresponds to A c = 0.498. 
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D. Infinite-State Truncation: Free Scalar Field Theory 



Lastly, we return to the case of free scalar field theory in 1 + 1 dimensions, but this time, 
we use a truncation algorithm which keeps an infinite number of states at each step. 

Consider a truncation procedure based upon two-site blocks. For each two-site block B p , 
we introduce the operators 



$(p) ± = -^[0(2p)±0(2p+l)], 
n(p) ± = ^=[U(2p)±U(2p + l)], 

and define ladder operators A+{p) and A-(p) using 



(65) 
(66) 



$(p)± 
Q(p) ± 



1 




A(p)l + A(p) 
A{p)i - A{p) ' 



± 



(67) 
(68) 



where 7- = \/n 2 + 3 and 7+ = \J /j? + 1. In terms of these variables, the two-site Hamilto 
nian is simply a sum of two decoupled oscillators, and its eigenstates are given by 

1 



\N + (p),N.(p)) 



At(p) N +AHp) N -h + n-), 



(69) 



where ^4 + |7 + , 7_) = A-\^ + , 7_) = 0. We now adopt a simple truncation procedure in which 
we keep an infinite set of block-states 



\N(p)) = \N + (p)=N; N.(p) = 0). 



(70) 



In other words, only states for which the higher-frequency 7_ oscillator is in its ground state 
are retained. With this choice of eigenstates, hi(p) conn can be written as 



2 2 



A{p)A + {p) + \ 

n + (p) 2 + 1 l$ + (pf 



(71) 
(72) 



Now consider an r-block sublattice B(p, ■ ■ ■ , p + r — 1) = {2p, 2p + 1, • • • , 2p + 2r — 1}. 
The Hamiltonian restricted to this sublattice has the form 
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j 2p+2r-l i 2p+2r-l 

^W = 2 E n 0') 2 +o E *0'i) M-O'i.Ja) tth), (73) 
where M r is a 2r x 2r real-symmetric matrix whose elements satisfy M r (i,j) = Ai r (i — j). 
M r can be diagonalized to obtain the normal modes, and the r-lowest eigenvalues of M r 
then yield Hr, the remnant eigenvalues. The triangulation matrix S is determined as usual, 
except that we can now work in terms of the fields instead of basis states. The proof of these 
statements is a straightforward exercise in normal-ordering using simple generalizations of 
the identities given in Appendix |A] and the definition of the operator S. The connected 
range-r term h r (p) conn is then computed from H r (p) by subtracting from it the previously 
computed connected range- j terms, for j = 1, ... ,r — 1. Finally, the terms h r (p) conn are 
combined to form the renormalized Hamiltonian, which takes the form 



rrren \ * ^_ 

~ ^ 2 



r-l 

(74) 



h + Q + {p) 2 + off $+(p)$+(P + ™ 

m=0 

where ho and a£j are c-numbers. For example, if we truncate the cluster expansion after 
two-block clusters (r = 2), we find in each RG step n — > n + 1 that a^\n + l) = a^\n)/2 
and ac(?\n + l) = a^\n) + (v^5 — l)a ( i\n)/2, for oti\n) < 0. Amazingly, no $ 4 terms 
appear in the renormalized Hamiltonian, but note that the exact CORE transformation of 
the nearest-neighbor Hamiltonian results in a new Hamiltonian which has an infinite number 
of terms. The importance of these results is that it shows we can, both in principle and in 
practice, directly deal with field theories having an infinite number of states per site, without 
first mapping them to spin systems. A more detailed description of the above calculation 
will appear in a forthcoming paper, which will also describe the analogous calculation for 
Fermi fields. 

At this point, an interesting question is "How big must r be in order to do a good job 
of reproducing the mass gap and correlation functions of the free-field theory?" . To analyze 
this question for the massless field (the hardest case), expand the fields $ + and Q + in terms 
of their Fourier components to rewrite the renormalized Hamiltonian for \x = as 



" /2 1 /- ~ 

k=-n/2 ^ V 



r-l 



cos(sk) 

.s=0 



$(fc)4>(-jfe) . (75) 
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We then explicitly compute the couplings for various values of r, the cluster expansion 
truncation order. Table |1| compares the values of Oq \ ct(\ a^"\ and ajj obtained from a 
range-2, range-3, and range-4 CORE computation. We see from the table that any given 
coefficient converges rapidly in r to its r —>■ oo limit. In this limit, the renormalized Hamilto- 
nian matches the original theory restricted to the subspace spanned by the oscillators having 
momenta —tt/2 < k < 7r/2; hence, we can use 



^ cos(sk) = 1 - cos(A;/2), (76) 
to determine the couplings in the r — > oo limit. We find 

at } = 1 - -, (77) 

7T 

^'-^&y (•><» < 78 > 

Note that the exact coefficients a^ 00 -* fall off as 1/s 2 . This means that if we truncate the 
formula for the frequency 

r-l 

uf(k) =2J2 cos(ks), (79) 

then the mass m(r), defined by m 2 (r) = uj^(0), fails to vanish; in fact, for some values of 
r, it becomes negative; the behavior of the g function of r is shown in Fig.||]. Since 

negative values for the gap make no sense as they imply that the renormalized Hamiltonian 
has no ground state, we can only truncate after an even number of terms. An important 
observation we can make from the plot is that we can only accurately compute the mass 
gap for the free field theory if we work out to a range r pa 1/m. We should note, however, 
that the mass is the quantity which is most sensitive to making a finite-range truncation of 
the exact renormalized Hamiltonian and that the function u r (k) converges for k > m more 
rapidly than it does for k = 0. This means that even for the worst case of a free-field theory, 
careful examination of the behavior of the function u r (k) and the oscillatory behavior of 
m(r) will allow us extract the correct physics without having to compute an infinite number 
of terms in the finite-range cluster expansion. A final point which we state without proof is 
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that the importance of the longer range terms is significantly reduced if we work with larger 
blocks. In one dimension for blocks of size L, the corresponding coefficients a s (r) fall off as 



This discussion applies to the truncation procedure in which we keep an infinite set of 
states at each truncation step. It would be interesting to discover how much of this behavior 
occurs if we keep only a finite number of states at the first step and map the theory into a 
spin system. The results of the Ising calculation show that we do not always need to work 
with large values of r in order to correctly reproduce the mass gap for a theory near its 
critical point. 

Note that CORE's ability to reproduce the mass gap and density of states near zero mo- 
mentum is much greater than that of the naive (t = 0) renormalization group procedure. For 
example, if we diagonalize the single-site Hamiltonian of the massless free scalar field theory 
and keep the single-site ground state and first excited state, then the naively-determined 
renormalized Hamiltonian is the simple Ising Hamiltonian for a value of the coupling far 
from its critical point; the system has a nonvanishing mass gap. However, if we keep the 
same single-site states but use the CORE rules to construct the renormalized Hamiltonian, 
many more interaction terms emerge and the couplings in the renormalized Hamiltonian 
are such that the system is much closer to the critical point where the mass vanishes. A 
more complete analysis of this system which examines the costs and benefits of keeping more 
states versus computing longer-range connected contributions to the renormalized Hamilto- 
nian would be very informative. 



The definition of the renormalized Hamiltonian in Eq. [1] is the cornerstone of the CORE 
approach. We were led to this definition by first observing that the state 



l/r 2 L 2 . 



IV. DERIVING THE BASIC RULES 



e 



tH/2 





(80) 
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contracts, as t — > oo, onto the lowest energy eigenstate of H for which the starting trial 
state |^o) has a nonvanishing overlap, typically the ground state. Note that the ground- 
state expectation value of any operator O can be obtained by taking the limit 

(0) = lim(*(t)\Omt)). (81) 

The use of Eq. |T] to define the renormalized Hamiltonian is strongly suggested by the following 
theorem. 

Theorem I: Given a Hilbert space H and truncation algorithm with associated 
projection operator P, diagonalizing the renormalized Hamiltonian 

H^{t) = [T(t) 2 H [T(t)HT(t)] mtftK (82) 

where T(t) = e~ tH and [O] = P^OP, is equivalent to finding the vector in 
the subspace V = P7i which minimizes the ratio 

E^t) = (^\T(t)HT(tM)/^\T(t) 2 \^). (83) 

Proof: Let \<p n ) label a basis for V = PTi and expand \ip) = J2 n a n\4>n)- Then, 

\T(t) HT(t)\<f> n ) 

To minimize this expression over all the states in V, differentiate Eq. with 
respect to a* and equate to zero; this yields 

(<f> r \T{t)HT{t)\i/>) =S4t)(<f) r \T(t) 2 \^). (85) 

Since this applies for each r, this can be rewritten 

[T{t)HT{t)m=£^t){T{tfm. (86) 



Hence, finding the state which minimizes S^(t) is equivalent to solving a gen- 
eralized eigenvalue problem. Given that {T(t) 2 ] is a positive matrix, the relative 
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eigenvalue problem can be converted into an equivalent standard eigenvalue prob- 
lem. In other words, finding the state which minimizes the ground-state energy 
of the original Hamiltonian is equivalent to diagonalizing the operator H ren (t) 
defined in Eq. 

In order to demonstrate that the renormalized Hamiltonian given in Eq. [I] defines a 
valid renormalization group transformation H Ten = t(H), we must at least show that the 
eigenvalues of H mn and the low-lying eigenvalues of the original Hamiltonian H are the 
same. 

Theorem II: Let if be a Hamiltonian defined in an iV-dimensional Hilbert space 
Ti with eigenstates {|e,g)} and corresponding eigenvalues e@, for (5 = . . . N — 1 
and €p < e/3+i. Let V be an M-dimensional subspace V C TC spanned by the 
states for a = . . . M— 1, where M < N. The projector into this subspace 

is P = Y^a^o 1 Furthermore, let S denote the M x M unitary triangula- 
tion matrix whose construction has been previously described, and denote 
the remnant eigenstates of H in V with corresponding eigenvalues Tp. Then the 
operator defined by 

H rcn = lim [T(t) 2 ]- 1/2 \T(t) H T(t)i [T(t) 2 ]- 1/2 , (87) 

t — ^oo 

where the contractor T(t) = e~ tH and \0\ = P^OP, simplifies to 

H Ten = S^H r S, (88) 

where H T = diag(7^, . . . , T M -i). 
Proof: Define the M states |£ a ) by 

M-l 

= E \*p) $U- ( 89 ) 

From the construction of S and the remnant eigenstates using the singular value 
decomposition, the states |£ a ) satisfy (Tp\£ a ) = for p < a. However, since 
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the projections of missed or non-remnant eigenstates into V can be expressed in 
terms of the projections of lower-lying remnant states, this means that (e p \C, a ) = 
for all p satisfying e p < T a . The use of Eq. |b| and the singular value decompo- 
sition to define the remnant eigenstates also ensures that (T p \£ a ) = even for 
p corresponding to remnant states which are degenerate with but orthogonal to 
\T a ), and hence, 

e- tH \U = e- tT "\T a )(T a \Z a ) + £ \e„) (e^ a ) , (90) 
where (T a \£ a ) ^ 0. Now define a new set of states 

\e a (t)) = z aP {t) (91) 

Z^(t) = 5 aP e tr -\(Wr a )\- 1 . (92) 

It is not difficult to show that 

\im(Q a (t)\T(t) 2 \Qp(t)) = 5 a p, (93) 

t — >CO 

\im(e a (t)\T(t)HT(t)\e p (t)) = TJ a/3 . (94) 
In matrix notation, these equations are 



lim Z(t)S\T(tf\S ] Z(t) = I, (95) 
lim Z{t)S[T{t)HT{t)]S^Z{t) = H T , (96) 



t—rOO 



where I is the M x M identity matrix. It then follows from Eq. |95| that 

_ i 

Iimf5[T(t) 2 ]5t)" 3 = limZ(t), (97) 

t — >oq V / t— *-oo 



and thus, Eq. [96| becomes 



lim (S[T(t) 2 ]St ) ( S[T(t)HT(t)]& ) ( S[T(t) 2 ]St ) = Hr . (98) 
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Using the matrix relation BA X I 2 B 1 = [BAB x ) x l 2 and the unitarity of S, it 
then follows that 



hm SlT(t) 2 j-HT(t)HT(t)jlT(t)T h S^ = H T , (99) 

and finally, 

H ren = giff^ ( 100 ) 

This theorem demonstrates that the eigenvalues of the renormalized Hamiltonian are the 
M eigenvalues T a associated with the remnant eigenstates of the original Hamiltonian. If 
the truncation procedure is such that no eigenvalues are missed, then the eigenvalues of H rcn 
are the M lowest eigenvalues of H. By showing that the mapping r{H) = H ren replaces 
the original theory with a theory containing the same low-energy physics but defined in 
terms of fewer degrees of freedom, Theorem II provides the justification for identifying r as 
a renormalization group transformation. 

An important aspect of the CORE approach is the use of the finite cluster method 
to approximate the renormalized Hamiltonian on the infinite lattice. In Ref. lattice 



constant theory was used to show that the finite cluster method can be applied in the 
calculation of any quantity so long as that quantity is extensive. We now demonstrate the 
extensivity of the renormalized Hamiltonian. Recall that a quantity is extensive if, when 
evaluated on a disconnected graph, it is the sum of that quantity evaluated separately on 
the connected components of the graph. 

Theorem III: The renormalized Hamiltonian is extensive. 

Proof: Consider a disconnected sublattice G = G\ U G 2 comprised of two con- 
nected components G\ and G 2 . Since G is disconnected, H{G) = H{Gi) + H{G 2 ) 
and [H(G 1 ),H(G 2 )} = 0; hence, T G (t) = T Gl (t)T G2 (t) = T G2 (t)T Gl (t) . Since the 
truncation is done on a block-by-block basis, then [T G (t)J = [T Gl (t)] [T Ga ( 
and 
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H\G) = [T Gl (t) 2 ] _1 / 2 |T G2 (t) 2 ] -1 / 2 ^[T Gl {t)H{Gi)T Gl {t)\\T G2 {t) 2 \ 
+ [T G2 (t)/J(G 2 )T G2 (t)][T Gl (t) 2 ])[T Gl (t) 2 ]- 1 / 2 [T G2 (t) 2 ]- 1 / 2 , 
= [T Gl W 2 r 1/2 lT Gl (t)ff (Gi)T Gl (0] [T Gl (t) 2 ]- 1 / 2 
+ [T G2 (t) 2 ]- 1/2 [T G2 (t)ff(G 2 )T G2 (t)][T G2 (t) 2 ]- 1 / 2 , 



(101) 



(102) 



H t {G 1 ) + H t {G 2 ). 



(103) 



Hence, if r is extensive. 

Clearly, the eigenvalues of the renormalized Hamiltonian on a given cluster containing 
R blocks are the same as the lowest M R eigenvalues (modulo missing ones) of the full 
Hamiltonian restricted to the cluster, assuming M states are retained in each block. In 
truncating the cluster expansion of the renormalized Hamiltonian on the infinite lattice, 
the correspondence between the low-lying eigenvalues of the infinite lattice H Ten and H can 
then be only approximate. However, our previous examples suggest that truncating the 
cluster expansion of the renormalized Hamiltonian after only a very few terms can lead to 
remarkably accurate results. 



In this section, we discuss two issues related to approximations and the CORE procedure. 
First, within the context of free scalar field theory, we link the methods presented in this 
paper to our earlier Physical Review letter where we used approximate contractors to carry 
out the computations. We demonstrate how such approximations converge and show why 
previous approaches always found a best finite value of t for determining the ground-state 
energy. Second, we reconsider the question of single-state truncations in the Heisenberg 
antiferromagnet. We do this to show how simple single-state truncations can encounter 
problems with surface effects and how the multistate renormalization group algorithm avoids 
these problems. 



V. APPROXIMATION ISSUES 
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A. Connection To Earlier Methods: Approximate Contractors 



Our earlier version of the CORE procedure [jl|] used an approximate contractor T n (t) 
obtained by decomposing H = Hi + H 2 into two or more parts and writing 

T n {t) = [ e ~ tH ^ 2n e~ tH ^ n e- tH ^' 2n } " , (104) 

where Hi and H 2 are chosen such that e~ tHl and e~ tH2 could either be computed exactly 
or numerically to any desired degree of accuracy. The validity of this approximation follows 
from the fact that, for operators A and B, one can show that 

e S(A + B) = e 5A/2 e 5B/2 e C(S) £ SA/2 ? (1Q5) 

oo 

C(5) = J25 2j+1 2j+ i. (106) 

3=1 

In particular, for 5 = t/n we see that as n —>■ oo, the sequence T n (t) converges to e~ tH as 
(t/n) 3 . When approximating the contractor in this way, we will see that t must be viewed 
as a variational parameter to be optimized. We will also see that this earlier procedure is 
less accurate and more time consuming than the method presently proposed. 

To see this, consider once again the Hamiltonian given in Eq. [19| which describes a free 
scalar field theory in one spatial dimension. This Hamiltonian can be expressed as a sum of 
single-site operators and nearest-neighbor interactions: H = H + V, where 



tfo = E^[n(j) 2 + 7o 2 0(j) 2 ], (107) 

3 

f = -£M>0' + i)> (los) 



3 

and 70 = vV 2 + 2. The ground state of H is then a product of uncorrelated Gaussians: 



<0|v]> o > = n<0l7o(j)> = ne" 70 ^ )2/2 . (109) 



o -7o<A(i) 2 /2 

3 3 



Our aim is now two-fold: to demonstrate how to recover the ground state of the full theory 
by applying e~ tH to j^o) an d taking the limit t — > oo; and to determine how well T n (t) = 
Ti{t/n) n = { e - tH o/2n e -tv/n e -tH /2n-^n a pp rox j m ates e~ tH for finite values of n and t. 
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First, evaluate e tH \^o)- Introducing the Fourier transforms 4>(k) and tl(k), we obtain 

H = Y, H{k) = Y.\ [u(-k)u(k) + u 2 (k) 4>{-k)4>{k)] , (110) 



i /,• 2 



111^ 



H = £ ff (Ar) = £ - [flC-^flCAr) + 7o 2 , 
k k 1 

V = Y J V(k) = - £ cos(A;) #(-*)#(*), (112) 

|^o) = Il ex P [-7o0(-A;)0(A;)/2] , (113) 

k 



where u(k) = J /i 2 + 4 sin 2 (A;/2). Since (f>(—h) = (j>(k)^ and fl( — k) = fl(ky, we can decom- 
pose the fields in terms of their real and imaginary parts and restrict all sums to k > 0, 
handling the case k = separately. Since the H(k) mutually commute, then e~ |^o) can 
be written as a product over states labelled by the momentum k so that we can limit our 
attention to a single fc-mode without loss of generality Let \j(k)) denote a simple harmonic 
oscillator ground state of frequency j(k). Now apply Theorem A2 proven in Appendix 

e- tH ^\ 1 (k)) = A(k,t)\ 1 (k,t)), (114) 

where A(k,t) is a normalization factor and \ j(k, t)) is a simple harmonic oscillator state of 
frequency 

T (M)=«(*) ^ - T - w — — . (115) 

Thus, as t — > oo, the frequencies j(k,t) — > which means that the state e~ tH \^ ) 

converges (up to the normalization factor YikA(k,t)) to the true ground state of the lattice 
free field theory. Since the normalization factor cancels out in ratios such as 

(%\e-'»"He-'»l'^„) 1„ / 

(« |e-|* 0> = 2gr*''> + ^J- (U6) 

we will ignore it from here on. 

The determination of T n (t)\ty ) proceeds similarly to that of e~ tH \^o). Since we have 
[H (k),H (k')] = 0, [V(k),V(k')} = and [H (k), V{k')] = for k ^ k', then T n (t)\V ) can 
be written as a product over states labelled by the variable k and we can study the general 
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problem one fc-mode at a time. If \j p (k)) is a simple harmonic oscillator ground state of 
frequency j p (k) associated with a mode k, then 



e- tH ^' 2n \ lp {k)) oc | 7 '(M/2n)>, (117) 
e-' v '( fc )/"|7 P (A;))oc|7 /, (A:,t/n)), (118) 



where the frequencies of the new oscillator ground states are related to j p (k) by 

7 (fc, t = 70 7T\ ; _ 2trynf / 1 \ 7 , 119 

V7 P (^)+7o-e 2t ^{j p {k) -jo) J 
j"(k,t) = j p (k) -2tcos(fc). (120) 

Using the above relations, [e~ tHo ^ k >' 2n e~ tv ^ k >' n e~ tH °^ k '' 2n ] n \j (k)) can then be easily evalu- 
ated. We find that the Gaussian state |tq(A;)) evolves to a new Gaussian state \j n (k,t)) of 
frequency j n (k,t). 

Plots of j n (k, t) for various values of k and n and a range of t values are shown in Figs. |6|, 
[7| and |8|. Plots of the expectation value of H in the state T n (t)\ty Q ) for the same values of n 
and range of t are shown in Fig. There are two things to notice about these figures. First, 
for larger values of k, the frequencies converge quickly to the values they would have in the 
exact wavefunction, indicated by the horizontal lines; however, for smaller values of k, the 
exact frequencies are not well reproduced, even for very large values of n. This means that 
computing the action of T n (t) on |^ ) can do well at approximating the ground-state energy 
density and still fail to reproduce the mass gap. Second, we observe that for finite values of 
n, there is a finite t which yields a best estimate of the ground-state energy density. T n (t) 
does a very good job of approximating e~ tH for smaller values of t, so at first, T n (t)\^fo) tends 
towards the ground state of H; however, T n (t)\^?o) eventually begins to move away from the 
ground-state wavefunction and so the expectation value of the energy density starts to get 
worse. This shows that without additional improvements, working with T n (t) for finite n 
and a best t cannot be expected to always accurately reconstruct the infrared properties 
of the theory. The renormalization group method works better than simply evaluating the 
action of T n (t) on a single state because it eliminates only the higher states which T n (t) 
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reproduces well and carries the more difficult long- wavelength modes over to the next step 
of the calculation. The agreement between the results of our earlier CORE treatment of the 
Ising model which used n < 16 and a best value of t and our current n = oo and t = oo 
calculation supports this picture. 

B. Antiferromagnet: Simple Cluster Formulae 

We now return to the Heisenberg antiferromagnet and compute the vacuum energy den- 
sity using two different single-state truncation procedures. There are two reasons for doing 
this: first, to show that computing the ground-state energy density for an infinite-volume 
theory from a series of finite-volume calculations is generally applicable; second, we often 
learn more from examples which do not work as expected than from ones which work well. 
In this case, we will learn that partitioning the lattice into either two- or three-site blocks 
can produce sequences of truncated cluster expansions which converge at very different rates. 
We explain why this happens and show how the two-state truncation algorithm used earlier 
avoids these convergence problems. 

First, we apply a single-state RG algorithm in which the lattice is partitioned into two- 
site blocks and we retain only the lowest-lying eigenstate in each block. Denote by E r 
the ground-state energy of the theory defined by restricting the full Hamiltonian to an r- 
site sublattice. The two-, four-, six-, and eight-site ground-state energies are E 2 = —0.75, 
E 4 = -1.616025, E 6 = -2.493577, and E 8 = -3.374932, and they yield the following 
connected contributions in the cluster expansion of the renormalized Hamiltonian: 



62 


— E2 


= -0.75, 


(121) 


e 4 


= E 4 


-2e 2 = -0.116025, 


(122) 


ee 


= Eq 


-2e 4 -3e 2 = -0.011527, 


(123) 


es 


= E 8 


- 2e 6 - 3e 4 - 4e 2 = -0.003803. 


(124) 



Thus, we obtain a sequence of approximations to the infinite-volume ground-state energy 
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density from the following truncated cluster expansions: 

S 2 = e 2 /2 = -0.375, (125) 

£ 24 = ( e2 + e 4 )/2 = -0.4330125, (126) 

£ 246 = (e 2 + e 4 + e 6 )/2 = -0.438776, (127) 

^2468 = (e 2 + Q + e 6 + e 8 )/2 = -0.4406775, (128) 

which are to be compared to the exact energy density £ ex act — —0.443147. Note that we 



divide by two in the above formulas so that our results refer to the energy per site of the 
original lattice instead of the energy per two-site block. For this simple truncation algorithm, 
the finite-range cluster expansion converges rapidly and agreement with the exact answer to 
better than one percent is obtained with ease. Given our earlier discussion of the free-field 
theory, it is interesting to compare the approximations built from connected terms to what 
we would obtain from simply dividing the ground-state energy for each n-site block by n. 
The comparison of these results is presented in Table |TTT| . 

The better than one-percent agreement of the finite-range cluster expansion with the ex- 
act ground-state energy density brings into question the benefits of using the renormalization 
group algorithm. However, the need for the renormalization group becomes apparent after 
examining the sequence of approximations obtained using three-site blocks. In this case, 
numerical diagonalization of the appropriate sublattice Hamiltonians yields E3 = —1.0, 
E 6 = —2.493577, and Eg = —3.736322, yielding connected contributions 

e 3 = E 3 = -1.0, (129) 
e 6 = E 6 - 2e 3 = -0.493577, (130) 
e 9 = E 9 - 2e 6 - 3e 3 = 0.250832. (131) 

If we now use these results to construct the corresponding approximations to the energy 
density per site, we obtain the sequence 

£3 = e 3 /3 = -1/3, (132) 
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£ 36 = (eg + e 6 )/3 = -0.497859, (133) 
^369 = (e 3 + ee + e 9 )/3 = -0.4142483, (134) 

which oscillates about the correct answer and converges much more slowly than that for 
the two-site decomposition of H. The cause of this oscillation and slow convergence arises 
from the fact that the physical excitations of this model have integer spin; the three-site 
decomposition has difficulty reproducing the low-lying physics since the ground state of 
the three-site block is a spin-1/2 multiplet, that of the six-site block is spin-0, and the 
ground-state of the nine-site block is once again spin-1/2. The two-site decomposition 
of the Hamiltonian does not suffer from this effect. This lack of rapid convergence is very 
instructive; since there is no way to know in advance what the correct spectrum of excitations 
is, this shows that we need a method for summing, at least partially, an infinite number of 
terms in the finite-range cluster expansion. As we saw in our earlier discussion of the 
antiferromagnet, this is what the full renormalization group calculation allows us to do. 



VI. LOOKING AHEAD 

This paper sets forth the basic rules for CORE computations, derives the rules from first 
principles, and discusses issues related to the convergence of the procedure. Future papers 
will focus on the application of these methods to more interesting physical systems and on 
clarifying the connection of the CORE approach to perturbative methods in instances where 
both are applicable. Some systems which should receive early attention are lattice gauge 



theories with and without fermions, t-J models ||18|| , and extended Hubbard models fL9fl . It 
is important to study the application of CORE technology to lattice gauge theories in order 
to see if, as we believe, it provides a powerful alternative to Monte Carlo calculations for 
studying QCD and chiral symmetry breaking. Extended Hubbard and t-J models are of 
interest because they are conjectured to have some relevance to high-T c superconductivity 
and have proven difficult to study in more than one spatial dimension by conventional 
methods. In this section, we discuss the application of CORE methods to these problems 
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and indicate how one could establish the connection between the CORE approach and a 
perturbative renormalization group treatment of 4 theory. 

A. Lattice Gauge Theory Without Fermions 

There are many ways to apply the techniques introduced in this paper to lattice gauge 
theories. One interesting approach is to divide the lattice into finite-size blocks, truncate the 
Hilbert space associated with each block to a set of gauge-invariant states, and then use the 
renormalization group formalism to map the gauge theory into a system which, like a spin 
system, has only a finite number of states associated with each lattice site. This approach 
yields an "equivalent" Hamiltonian theory in which all of the unphysical degrees of freedom 
have been eliminated. We can then treat the new Hamiltonian in the same way as in the 
Heisenberg and Ising models. 

For example, we could associate with each plaquette of the original lattice a single site 
in the new lattice. We could then find the low-lying gauge-invariant eigenstates of the 
one-plaquette Hamiltonian, either exactly or numerically, and truncate by selecting a finite 
number of these eigenstates. Using this truncation procedure, we construct a renormaliza- 
tion group transformation which maps the gauge theory into a generalized "spin" system. 
The interactions between nearby "spins" are found by evaluating the renormalized Hamil- 
tonian on clusters containing several connected plaquettes. This new spin system would be 
guaranteed to have the same low-lying gauge-invariant physics as the original theory and 
could be treated in the same way as the Heisenberg and Ising models. This approach allows 
us to define and carry out a gauge-invariant renormalization group calculation for any lattice 
gauge theory. 

This ability to define a gauge- invariant, Hamiltonian-based, real-space renormalization 
group calculation is unique to the CORE approach. Earlier real-space renormalization group 
procedures also kept a finite number of states per block, but they defined the renormalized 
Hamiltonian by [[H]], the truncation of the original Hamiltonian to the subspace spanned 
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by the retained states (this corresponds to the t = limit of the CORE approach). In such 
calculations, keeping only gauge-invariant block states leads to a truncated Hamiltonian 
in which the block-block interactions vanish. In order to retain inter-block couplings, flux 
must move across the links joining the blocks; this cannot happen without keeping some 
gauge-noninvariant single-block states. However, if one keeps such states in the truncation 
procedure, the entire process becomes much more cumbersome. 

The question of how many single-block gauge-invariant states and how many terms in 
the cluster expansion of the renormalized Hamiltonian should be retained naturally arises 
when carrying out a contractor renormalization group calculation; each choice constructs a 
mapping of the original gauge theory into a different generalized spin system. We hope to 
answer this question in the future by carrying out several computations in a simple lattice 
gauge theory, such as 2+1- dimensional compact U(l), varying the number of retained single- 
block states and clusters to see how quantities of interest, such as mass gaps and the specific 
heat, depend on these factors. 



B. Lattice Gauge Theory With Fermions 

Interesting possibilities arise when we consider lattice gauge theories with fermions. One 
way of treating these theories is to study systems with either SLAC [pO[j , Wilson or 



Quinn-Weinstein []22| fermions and truncate the system to the subspace spanned by tensor 
products of gauge- invariant, single-site states. In the case of lattice QCD, this would in- 
clude all color-singlet single-site states, i.e., mesons and baryons, which can be formed by 
applying quark and antiquark creation operators to the single-site vacuum state, subject to 
the constraints imposed by the exclusion principle. As the only terms which appear in the 
lattice QCD Hamiltonian create (or destroy) closed loops of flux or move quarks from site 
to site trailing their flux behind them, the color-singlet mesons and baryons are all degen- 
erate and the connected range-1 part of the renormalized QCD Hamiltonian will vanish. In 
order to compute the connected range-2 terms, we solve the problem of two sites connected 
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by a single link and find the low- lying gauge-invariant eigenstates which have an overlap 
with all of the tensor products of the two sets of single-site meson and baryon states. This 
computation yields connected range-2 contributions to the renormalized Hamiltonian which 
contain meson and baryon kinetic terms as well as meson-meson and meson-baryon interac- 
tions. Connected range-3 terms come from computations involving three sites arranged in a 
straight line or forming a right angle. These range-3 terms contain corrections to the terms 
already described, new terms which allow mesons and baryons to hop along diagonals of 
the underlying lattice, and terms which describe three-site interactions. Continuing in this 
way produces a renormalized Hamiltonian expressed only in terms of the physical degrees 
of freedom; the underlying quarks and gluons disappear from the problem. 

We would now like to say something about how chiral symmetry breaking will show up in 
QCD with three flavors of quarks and either SLAC or Quinn-Weinstein fermion derivatives 
(the case of Wilson fermions is somewhat different). Consider a theory with three massless 
flavors of quarks and apply a more restrictive truncation procedure which keeps only single- 
site fluxless states containing equal numbers of quarks and antiquarks, i.e., mesons. For 
three flavors of quarks there are 924 such states and, as was shown in Ref. [33[], they form an 



irreducible representation of the group SU(12) where the group generators are formed from 
bilinears in the single-site quark fields Q. Note that for Nf flavors, the fluxless states form 
an irreducible representation of the group SU(4Nf). For a truncation algorithm based upon 
keeping gauge-invariant single-site states, the renormalized Hamiltonian contains no range-1 
connected terms. The first nonvanishing contribution to the renormalized Hamiltonian will 
be the range-2 connected terms and these are computed by solving the two-site theory. 



It was pointed out in Ref. |23| that if we keep only the nearest-neighbor terms in the 
fermion derivative, then the resulting Hamiltonian is invariant under a global SU(12), and 
since the two-site problem cannot have anything but nearest-neighbor terms, this observation 
can be used to simplify the computation of the connected range-2 terms in the renormalized 
Hamiltonian. We already noted that the fluxless single-site states form an irreducible 924- 
dimensional representation of SU(12) and so tensor products formed from these states can 

44 



be decomposed into the irreducible representations of SU (12) which appear in the product 
of two 924's; these are the only states in the full problem relevant to our CORE compu- 
tation. Starting from the highest weight state in each of these irreducible representations 
and applying the Lanczos method, we can numerically find the relevant eigenvalues of the 
two-site Hamiltonian to a high degree of accuracy. From general symmetry arguments, the 
most general two-site Hamiltonian one can write for this system will be in the form of a finite 
polynomial in the Casimir operator and higher order invariants formed out of the generators 
of SU (12). Thus, the general structure of the connected range-2 Hamiltonian will be given 
by 

hT a U) = ocxQ(j) ■ QU + 1) + <*2(Q(j) ■ Q(J + I)) 2 + • • • • (135) 

It is a simple exercise to show that in the strong-coupling limit, the leading term in this 
expansion is the one proportional to Q(j) ■ Q(j + 1); in other words, in strong-coupling, the 
renormalized range-2 Hamiltonian is just a generalized Heisenberg antiferromagnet. As was 
argued in Ref. [^3] and Ref. |2"4| , we expect this theory to spontaneously break to SUy(6) x 



SUa(Q), where the vector SUy(6) is realized normally and the axial-vector SUa(Q) is realized 
in the Goldstone mode. Thus, in the strong-coupling limit, the connected range-2 part of the 
renormalized Hamiltonian unavoidably leads to a spontaneously broken symmetry, but the 
group is too large and there are too many Goldstone bosons. Clearly, a detailed calculation 
is necessary to determine if these conclusions persist in weak coupling where other terms 
in Eq. |135| can become significant. However, we can show that the problems of having too 



large a symmetry group and too many Goldstone bosons disappears once we compute the 
connected range-3 terms. 

To see this, observe that, independent of the coupling constant, the next-to-nearest- 
neighbor terms in both the SLAC and Quinn-Weinstein types of derivative break the SU (12) 
symmetry and, after including these terms in the renormalized Hamiltonian, all that remains 
of the SU(6) x SU(6) symmetry of the nearest-neighbor theory is SU(3) x SU(3). As in 
the discussion of the range-2 terms, we can invoke the strong-coupling limit to calculate the 
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structure of the leading range-3 terms and explicitly show that the range-3 terms give the 
unwanted Goldstone bosons mass and that the degenerate SU(6) multiplets of mesons break 
up into SU(3) multiplets. This is in strict analogy to what was discussed in Ref. [p4"| . Of 
course, as we noted for the case of the range-2 terms, the generic structure of the connected 
range-3 terms in the renormalized Hamiltonian is richer than that of the leading terms in 
the strong coupling limit, and so asserting that this pattern of symmetry breaking persists 
to the physically more interesting weak-coupling regime requires more work than we have 
done to this point. 

Much interesting work remains to be done in this picture of dynamical chiral symmetry 
breaking; nevertheless, the fact that the CORE procedure provides a coupling-independent 
way of constructing an effective theory of mesons which, in the strong-coupling limit, coin- 
cides with earlier descriptions in which dynamical chiral symmetry breaking appears natu- 
rally, is new and unique to this approach. 

C. Hubbard and Extended Hubbard Models 

Among the interesting features of the Hubbard and extended Hubbard models are the 
variety of phase transitions which can occur as the density of particles in the ground state 
changes. While tuning the density of particles in the ground state is easily accomplished 
by adding a chemical potential to the Hamiltonian, early attempts to analyze these theories 
using naive real-space renormalization group methods ran into problems: projecting onto 
a small number of states per block so that the occupation number of each state is a finite 
integer, and therefore the density a rational fraction, made it difficult to achieve a smooth 
dependence of the density on the chemical potential. CORE mitigates this problem without 
having to keep a large number of states per block: first, the connected range-r terms are 
computed by diagonalizing the full r-site Hamiltonian, including the chemical potential, and 
so these terms can encode more complicated behavior of the chemical potential coefficient 
fjL] second, the operator which measures the density of particles in the ground state as a 
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function of \i undergoes a much more complicated evolution than it does in a naive truncation 
procedure, evolving connected range-r terms of its own. Preliminary computations support 
this picture but more extensive computations are needed to fully explore the potential of 
CORE methods for this class of problems. 



D. Connection to Perturbation Theory 

In this section, we discuss the way in which one could establish the relationship be- 
tween the CORE approach and the familiar perturbative renormalization group in the weak- 
coupling limit. To illustrate this connection, consider adding a A0 4 interaction to the scalar 
field theory Hamiltonian given in Eq. [19], where the A coupling is small, and again apply the 
CORE procedure outlined in Sec. IlfTJ, It is a straightforward exercise to include the 4 
term and perturbatively compute the CORE transformation associated with the two-site, 
infinite-state truncation procedure. 

We begin with the same truncation procedure defined for the free-field case and keep 
the same tower of oscillator states for each two-site block. A new feature is that we must 
now compute S and Hq- even for the range-1 terms because the retained states contract 
onto states which are different from the two-site, free-field eigenstates and the eigenenergies 
corresponding to these states are also changed from their free-field values. A direct con- 
sequence of this is that the new range-1 connected part of the renormalized Hamiltonian 
contains higher-order polynomials in the fields. Given a perturbative expression for the 
connected range-1 terms, we have to perturbatively solve the four-site problem to compute 
S and H? in order to obtain the connected range-2 terms in the renormalized Hamiltonian. 
Once again, we get a set of terms of the form <p m (p)(f) n (p + 1) which do not correspond to 
terms in the original Hamiltonian. Longer-range connected terms are computed in the same 
manner. Since the zero-coupling limit of this procedure builds up a finite-range expansion 
of the free-field theory, one should be able to make this perturbative expansion match up 
with more familiar renormalization group computations. 
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There is a simple way to modify the procedure just outlined so as to automatically 
resum the perturbative expansion of the renormalized Hamiltonian to very high order in 
the coupling. One virtue of this modified approach is that it guarantees that the ground- 
state energy density will behave as A 1 / 3 for large couplings. The basic idea is to change the 
definitions of 7± in Eq. ^ in order to treat them as variational parameters which depend 
upon /i 2 and A. To determine their values, we minimize the expectation value of the two- 
site Hamiltonian in the state I7+, 7_) with respect to 7+ and 7_. Fixing 7+ and 7_ in 
this way, we then rewrite the two-site Hamiltonian in terms of annihilation and creation 
operators, normal order the resulting expression, and do perturbation theory in the non- 
quadratic terms. Note that this minimization process guarantees that the state |7 + ,7_) 
is the lowest-lying eigenstate of the "free Hamiltonian" obtained by keeping the quadratic 
terms, including those which come from normal ordering the quartic self-interaction. Since 
7+ and 7_ are nontrivial functions of fi 2 and A, the perturbation theory just described 
amounts to an infinite resummation of the usual expansion. 

To compute the range-2 terms in the renormalized Hamiltonian, solve the four-site free 
problem but treat the 7 frequencies as variational parameters determined by minimizing the 
expectation value of the Hamiltonian in the ground state of the oscillators. This leads to four 
coupled equations which can be solved numerically for any value of A. Once again, normal 
order the Hamiltonian and treat all terms which are not quadratic in the ladder operators 
as perturbations. The states obtained by working to finite order in these perturbations are 
used to construct S. The computation of higher-range connected terms proceeds in a similar 
manner. 

VII. CONCLUSION 

The contractor renormalization group, a general method for solving any Hamiltonian 
lattice system, was presented. The CORE approach is a systematic and nonperturbative 
procedure for carrying out real-space renormalization group transformations which relies on 
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contraction and cluster techniques. The method was illustrated using four examples: free 
scalar field theory with single-state truncation, the Heisenberg antiferromagnetic spin chain 
with two-state truncation, the anisotropic Ising model with two-state truncation, and free 
scalar field theory with an infinite-state truncation scheme. The use of approximate contrac- 
tors, the convergence of the cluster expansion in determining the renormalized Hamiltonian 
and the need for summation via the renormalization group were also discussed. 

A particularly exciting feature of the CORE technology is its ability to treat systems with 
dynamical fermions, systems which are difficult to study using stochastic methods. CORE 
also makes possible gauge-invariant renormalization group transformations in Hamiltonian 
lattice gauge theory and easily incorporates a chemical potential. These features suggest 
that the CORE approximation will prove to be a powerful tool in future applications to the 
Hubbard and t-J models and lattice gauge theory with and without fermions. 

This work was supported by the U. S. DOE, Contract No. DE-AC03-76SF00515 and 
Grant DE-FG03-90ER40546, and the UK PPARC through grant GR/J 21347. 
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APPENDIX A: HARMONIC OSCILLATOR IDENTITIES 

Consider the canonically conjugate operators x and p which satisfy [x,p\ = i, and intro- 
duce the one-parameter family of annihilation and creation operators A w = xJu/2+ip/ \/2uj 



and Al = xJuj/2-ip/V2u which satisfy [A*, A*] = 1- Note that x = (A~l + A^) / y/2uj and 



p = i (Al — A LU )Ju/2. Furthermore, define the Hamiltonian H u = \(p 2 + oj 1 x 2 ). 




Theorem Al: The state \u)\) defined by A^Auji) = and the state \lu ) defined 
by A^Auoq) = are related by 




(Al) 



Proof: Write A Wl in terms of x and p, then express x and p in terms of A Wo and 
A* to show that 



A m = 7m Al + 7p A 



(A2) 



where 




(A3) 



(A4) 



Now use [Au , A}?] = nA^ 1 to show that 




(A5) 



Hence, 
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exp 

(C4o) N> = (7m + 2^7p) exp (e4o) 4>o>, (A6) 

which vanishes if we set £ = — 7m/(27p) = (^ — c<j 1 )/[2(c<j +c<j 1 )]. Since A W1 = 
defines |c<Ji), then clearly 

N«i exp {-7M</(2 7 p)}ko) = ki). (A7) 
Requiring (cji|c<ji) = 1 and using {u \A™ A^™\u) ) = 5 mn n\, we have 



oo 



^ 2 = E^ 2n = (l-4a- 1/2 , (A8) 

n=0 V i V 

so = [4cjicj /(cj! +cj ) 2 ] 1/4 - 

Given this result, we can now easily show that applying e~ tHw ° to an arbitrary Gaussian 
wavefunction produces a new Gaussian wavefunction of a different frequency; as t — > oo, the 
new frequency tends to u . 

Theorem A2: Let |a^i) be the simple harmonic oscillator ground state defined 
by A^ui) = and H Uo be the Hamiltonian for a simple harmonic oscillator of 
frequency u . Then 



-tH, 



*\u 1 )=A(t)\u(t)), (A9) 



where \uj(t)) is the ground state of a simple harmonic oscillator of frequency u(t), 
and 

Proof- Using e~* Ht "(> A*™e tH "° = e~ mtu} ° A%£, one sees that 

e -*^o e £<e*^o = e ^o, (A12) 



53 



where £' = e 2fc<; o£. Using e |u; ) = e taJ °/ 2 |c^ ) and Eq. [Al], one finds that 

e-^o | Wl ) = Nui e-*»fl exp 1 e^A* } |u, >. (A13) 



If we set 



(cup - u>(t)) 
2(cu(t)+a;o) 



— 2 t lt)Q 



2 (wo + ^l) 



(A14) 



then we can identify the state on the right-hand side of Eq. [A13| with the Gaussian 
wavefunction \co(t)). Solving Eq. |A14| for u(t) yields the result given in Eq. |A11| . 
Note that uj(0) = uj\ and u(t — > oo) = u . The multiplicative factor A(t) is 
then given by A(t) = (N Ul / iV^Oe - *" / 2 which simplifies to the result shown in 
Eq. |AL0. 
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FIG. 1. Fractional error in CORE estimates of the ground-state energy density in the Ising 
model against A. The dotted curve with crosses shows previous estimates from Ref. results 
from the present work are shown by the solid curve with circles. 
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FIG. 2. CORE estimates (circles) of the mass gap in the Ising model against A. The solid curve 
shows the exact mass gap. 
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FIG. 3. Comparison of the CORE estimates (crosses) of the magnetization with the exact 
results (solid curve) in the Ising model against A. A c is the critical point. 
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FIG. 4. Plots of p(X) as given in Eq. 64 for A c = 0.496 (dashed curve with squares), 0.498 (solid 
curve with circles), and 0.500 (dotted curve with diamonds). A c is the critical coupling in the Ising 
model, and £ is the critical exponent corresponding to the magnetization. 
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FIG. 5. CORE estimates of the mass gap squared m 2 (r) in the free scalar field theory against 
the truncation order r in the cluster expansion of the renormalized Hamiltonian. The CORE 
estimates are obtained using an infinite-state truncation algorithm. 
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FIG. 6. The scalar field theory frequencies j n (k,t) obtained using an approximate contractor 
T n (t) for n = 1 and various momenta k. The starting state is a product of uncorrelated Gaus- 
sians. The curves correspond to different values of k; the corresponding frequencies in the exact 
wavefunction are indicated by the horizontal lines. 
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FIG. 7. The frequencies j n (k,t) as in Fig. ^, except that 
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FIG. 8. The frequencies 7 n (fc, i) as in Fig. ||, except that 



62 



FIG. 9. Expectation value of the free scalar field theory Hamiltonian in the state T n (t) |^o) for 
n = 1, 2, and 40, showing how the best t changes as a function of n. 
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TABLES 

TABLE I. Convergence of the partial sums e n in the cluster expansion of the free massless 
scalar-field vacuum energy density as a function of the range n. The energies are CORE estimates 
from a single-state per site truncation algorithm, and the errors are the differences of these estimates 
from the exact energy density 0.636619772. 



n 




Error 


1 


n 7A71 H7 
(J. fU71U7 


0.0 1 04o r 


Z 


U.oooyiy 


n nooono 

u.uzzzyy 


Q 




n £. A IP. A A 

U.04 / 044 


n ni 1 ni£ 
U.UllUzo 


A 

4 


U.t)4ozUo 


U.UUoooD 





U.u^lUUl 


U.UU400Z 


6 


0.639746 


0.003126 


7 


0.638962 


0.002343 


8 


0.638441 


0.001821 


9 


0.638076 


0.001456 


10 


0.637811 


0.001191 


20 


0.636932 


0.000312 


30 


0.636761 


0.000141 


40 


0.636700 


0.000080 
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(r) 

TABLE II. Comparison of some of the couplings a s in the renormalized Hamiltonian of a free 
scalar field theory for r = 2,3, 4, oo, where r is the cluster-expansion truncation order. Note that 
only ao and at\ are defined for the case r = 2. 



(r) 

as 



a 

Q 3 



r = 2 
0.381966 
-0.500000 



r = 3 
0.371054 
-0.451922 
0.104212 



r = 4 

0.367594 
-0.438360 

0.098137 
-0.041400 



r = oo 
0.363380 

-0.424413 
0.084883 

-0.036378 



TABLE III. Comparison of truncated-cluster £ r and finite- volume E r /r Heisenberg antiferro- 
magnet ground-state energy densities for range r. 



r £ r E r /r 

2 -0.3750000 -0.3750000 

4 -0.4330125 -0.4040063 

6 -0.4387760 -0.4155962 

8 -0.4406775 -0.4218665 
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